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ABSTRACT 

We investigate the effect of a distorted neutron star dipole magnetic field on 
pulsar pair cascade multiplicity and pair death lines. Using a simple model for a 
distorted dipole field that produces an offset polar cap, we derive the accelerating 
electric field above the polar cap in space charge limited flow. We find that 
even a modest azimuthally asymmetric distortion can significantly increase the 
accelerating electric field on one side of the polar cap and, combined with a 
smaller field line radius of curvature, leads to larger pair multiplicity. The death 
line for producing pairs by curvature radiation moves downward in the P-P 
diagram, allowing high pair multiplicities in a larger percentage of the radio 
pulsar population. These results could have important implications for the radio 
pulsar population, high energy pulsed emission and the pulsar contribution to 
cosmic ray positrons. 

Subject headings: pulsars: general — stars: neutron 

1. INTRODUCTION 

Rotation-powered pulsars are thought to produce electron-positron pairs through elec- 
tromagnetic cascades in their magnetospheres (Sturrock 1971). Such pair cascades may 
occur both in the magnetic polar regions, where high energy photons radiated by acceler- 
ated particles undergo pair conversion in the strong near-surface magnetic fields (Ruderman 
& Sutherland 1975, Daugherty & Harding 1982) or in outer gaps, where high energy photons 
interact with thermal X-rays from the neutron star (NS) surface (Cheng, Ho & Ruderman 
1986, Hirotani 2008). For some time though, there has been observational evidence that 
the number of pairs that can be generated in pair cascades by the existing standard models 
of particle acceleration in magnetic dipole fields is insufficient to account for the optical to 
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X-ray emission from the synchrotron nebulae powered by the pulsars. Estimates of the pair 
multiplicity (the number of pairs produced by each primary accelerated particle) of about 
10 5 — 10 6 needed to account for the emission from the Crab pulsar wind nebula (PWN) 
(DeJager et al. 1996) of about 10 5 for the Vela PWN (DeJager 2007) are more than one 
order of magnitude larger than the theoretical pair multiplicities. Synchrotron absorption 
models for the eclipse in the double pulsar system PSR J0737-3039 (Arons et al. 2005, Lyu- 
tikov 2004) require a pair multiplicity of around 10 6 for the recycled 22 ms pulsar in that 
system. There has also been a long-standing problem in understanding how long-period pul- 
sars are able to produce coherent radio emission, thought to require electron-positron pairs 
( Arons 1983). Again, the standard pulsar models are not able to account for the operation 
of robust pair cascades in these aging pulsars (Hibschmann & Arons 2001). More recently, 
the discovery of pulsed gamma-ray emission from a large number of millisecond pulsars by 
the Large Area Telescope (LAT) on the Fermi Gamma-Ray Space Telescope (Abdo et al. 
2009) has revealed light curves that are best modeled by narrow radiation gaps in the outer 
magnetosphere (Venter et al. 2009, Abdo et al. 2010). Such narrow gaps require screening 
of the accelerating electric field over most of the magnetosphere by a pair multiplicity that is 
orders of magnitude larger than standard models are able to produce (Harding & Muslimov 
2002). 

This motivates us to investigate the effect of a large-scale non-dipolar NS magnetic 
field geometry on pulsar pair cascades. Off-centered dipole-like magnetic fields seem to be 
prevalent among planets in the Solar System (Russell & Dougherty 2010), and in some 
stars. For example, the magnetic fields of isolated magnetic white dwarfs may possess an 
off-centered dipole or dipole plus quadrupole configurations (see e.g. Putney & Jordan 1995) 
which is consistent with the idea that the magnetic Ap and Bp stars, the predecessors of 
highly magnetic white dwarfs, favor an off-centered dipole field. However, it is very likely 
that the surface magnetic fields of white dwarfs are much more complex (Wickramasinghe 
2001). Recent modeling of X-ray pulse profiles from millisecond pulsars (MSPs) shows 
evidence for offset dipole fields or offset polar caps in PSR J0437-4715 (Bogdanov et al. 
2007) and PSR J0030+0451 (Bogdanov k Grindlay 2009). Modeling of the X-ray pulsations 
of the NSs in Low-Mass X-ray Binaries, the progenitors of millisecond pulsars, show possible 
evidence of even more extreme magnetic field distortions (Lamb et al. 2009), that could result 
from distortion of the global magnetic field by e.g. the crustal plate tectonics (Ruderman, 
1991). In pulsar theory, non-dipolar magnetic fields have been discussed by many authors 
and long before the Fermi era. Arons (1983) argued that pair creation in long period 
(P > 1 s) pulsars was difficult and that distortions of polar field lines due to large-scale, non- 
dipolar surface fields of magnitude comparable to the dipole component, and of much smaller 
radius of curvature, could enable pair creation in older pulsars. Arons (1997) concluded that 
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"while the frame dragging effect in star centered dipole geometry does improve comparison 
of the theory with observation, an unacceptably large fraction of the observed stars outside 
the bounds of pair creation theory still persists". He conjectured that dipole offsetting 
improves the correspondence between theory and observation, and could account for the 
"death valley" observed for pulsars approaching the radio death line. In this Letter we report 
our preliminary assessment of the effect of a distorted dipole magnetic field of a NS on the 
efficiency of pair creation. A more detailed treatment will appear in a later publication. 



2. Magnetic and Electric Field of a Distorted Dipole 

We propose the following heuristic model of an asymmetric magnetic field of a pulsar 
that can be used to approximate plausible distortions of the magnetic field of a NS. We 
introduce an azimuthal asymmetry to the field lines of an originally symmetric dipole such 
that the field lines over half of the polar cap (PC) have relatively smaller radius of curvature 
and over the other half of the PC have larger radius of curvature (see Figure 1 below). 
Consequently, one side of the PC is larger and the PC is effectively shifted from the [star] 
center of symmetry. 



In magnetic spherical polar coordinates (9, (ft) the magnetic field is assumed to have the 
following form: 

r„ r 1 . i 

(i) 



B ~ — — 

rj 6 



r cos[0(l + a)} + i 0sin[0(l + a)} 



where rj = r/R ns is the dimensionless radial coordinate in units of stellar radius, a is the 
parameter characterizing the distortion of polar field lines. Note that vector B is only 
approximately solenoidal, since we omitted a small 5^- component (~ B 9 da/d(f>) and higher- 
order corrections in B r and Bg that enter only higher-order terms neglected in the formulae 
below. In subsequent studies we will discuss these terms as well as more general types of 
field distortion. If the magnetic axis lies in the x-z plane, a = e cos produces an effective 
offset of the PC in the x-z plane and a = e sin <fi an effective offset of the PC in the y-z plane, 
where < e < 1. By using formula ([1]) we get the equation of the magnetic field line, which 
we give here in the small-angle approximation valid near the NS surface, 

9 = x (1+a ^ 2 , (2) 

1 + a 

where x = t/tlc is the radial distance in units of the light cylinder radius, tlc = c/Q; and 
< £ < 1 is the colatitude of a footpoint of polar field line normalized by the colatitude of 
the PC boundary. Figure 1 shows the field lines computed from Eqn (2) plotted in the x-z 
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plane. The field line radius of curvature (in units of tlc) is 
where 

a = sA, /3 = -A, A = eV +Q . (4) 
8 6 

The effective shift of the PC on the NS surface is approximately, 

Ar PC ~ R ns 9 [1 - 6 e Q ] , (5) 

where 9q = (fii? ns /c) 1//2 is the canonical half-angle of the PC, Q is the pulsar rotation rate 
and R ns is the NS radius. Since the effective offset of the PC is a fraction of the PC opening 
angle, it is a very small fraction of the stellar radius for normal pulsars and somewhat larger 
fraction for millisecond pulsars. In the rest of the paper we will be referring to e as the offset 
parameter. 

Neglecting the static general-relativistic corrections, we can write the accelerating elec- 
tric field, assuming the boundary conditions of space-charge limited flow (SCLF) (cf. Harding 
& Muslimov 1998), as 

E ""-2\—) Bo (rTWW [ +a{r] " 1)]cosx + 8(i-^)^ smxcos0 j (1 -°' 

(6) 

where x is the pulsar obliquity, k ~ 0.15 I45/RI (-^45 = -^/10 45 g-cm 2 , R 6 = i? ns /10 6 cm, I is 
NS moment of inertia) is the parameter accounting for the general-relativistic frame dragging 
( Muslimov & Tsygan 1992). Note that the second term in Eqn fl6]) is significantly simplified 
and is not valid for very large values of x- We also note that the current of primary charges 
assumed in the SCLF models, which is j = cpcj(O), is not consistent with the polar cap 
current that has been derived in force-free magnetosphere models (e.g. Contopoulos et al. 
1999, Timokhin 2006, Spitkovsky 2006), where Pgj(0) = QB /2nc is the Goldreich- Julian 
charge density at the surface. The discharge properties under these conditions were discussed 
by Beloborodov (2008) and may involve non-steady charge outflow and electric fields. We 
use here the SCLF electric field for comparison with the many previous studies of pulsar 
pair cascades and pair death lines. 



3. Pair Cascade Simulation and Results 



Using the distorted field structure and accelerating electric field described above, we have 
simulated a full pair cascade above the pulsar polar cap to compare the pair multiplicity 
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and pair death lines for different offset parameters e. The full pair cascade simulation 
code, described in Daugherty & Harding (1982) and Daugherty & Harding (1996), has been 
combined with the acceleration, early cascade development and pair front formation used in 
Harding & Muslimov (1998). The combined code is thus able to follow a primary particle 
through its acceleration in the parallel electric field, taking into account inverse Compton 
and curvature radiation losses, the emission of photons by curvature radiation (CR), the 
establishment of the pair front (above which the electric field is screened) and the full cascade 
of pairs and their synchrotron radiation. Since the electric field is screened over a very short 
length scale by CR pairs (Harding & Muslimov 2001), for the present calculation we simply 
assume that the electric field is zero above the pair formation front (PFF). Since the results 
are fairly sensitive to NS equation of state (EoS), we choose for normal pulsars a mass of 
Mq = 1.45 and radius R QS = 10 km, which give a moment of inertia / = 1.13 x 10 45 gcm 2 
(Lattimer & Prakash 2007), and for MSPs, that may have accreted mass during their spin- up 
phase, we choose M = 2.15, R ns = 9.9 km and / = 1.56 x 10 gem 2 from a rotating NS 
model (Friedman et al. 1986). 

The resulting pair multiplicities for a range of pulsar periods and surface magnetic field 
strengths, and for offset parameters e = 0, 0.1, 0.2, 0.4 and 0.6 are shown in Table 1 for the 
normal pulsar population and in Table 2 for millisecond pulsars. For these calculations, 
we have taken a = esin0 for an effective offset in the y-z plane and magnetic azimuth 
angle = 270°, where the field distortion produces the most favorable conditions for pair 
creation, i.e. smallest radius curvature, largest PC angle and largest electric field. Pair 
multiplicities are computed for one representative magnetic colatitude of £ = 0.5. The 
increase in the parallel electric field, in this case on the trailing side of the PC, by a factor 
of about 2a /(l + a) 2 , results in a larger voltage at the PFF and is the most important 
effect leading to larger pair multiplicities since the maximum curvature radiation energy is 
proportional to the cube of the particle Lorentz factor. For normal pulsars, even modest 
values of e that are equivalent to offsets of a few percent of the stellar radius lead to orders 
of magnitude increases in the pair multiplicity. Small offsets also enable long period pulsars, 
that cannot produce any pairs from CR in a pure dipole field, to produce high-multiplicity 
pair cascades. Most MSPs cannot produce pairs from CR in a dipole field but with some 
offset of the PC would be able to produce high pair multiplicity. The plateaus and peaks 
in the multiplicity, that can be seen in Table 1 for Bq > 3 x 10 12 G, are the result of pairs 
being produced increasingly in the ground Landau state in high fields (Baring & Harding 
2001) which inhibits the cascade synchrotron radiation. Photon splitting, which we have not 
included in these calculations, would further reduce pair multiplicity above 10 13 G. Although 
the multiplicity grows with increasing e and increasing field strength, it saturates at a value 
below 10 5 . This limit, also noted by Medin & Lai (2010), results from pairs being produced 
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with higher and higher average energy at smaller angles to the magnetic field as the primary 
electron energy, and thus the maximum energy of the CR photons, increases. 

The increase in pair multiplicity from a distorted dipole field can move many pulsars 
above the death line for robust pair creation. Figure 2 shows the death lines for pair creation 
by CR in the P-P diagram for different offset parameters, computed for the two different 
EoS for normal and MSPs described above. The lines plotted in Figure 2, determined by the 
method described in Harding & Muslimov (2002), show the lowest P required for a pulsar at 
a given period to establish a PFF. The death line for a pure dipole cuts through the middle 
of the normal pulsar population and through the upper edge of the MSP population. As 
e increases, the lines move significantly downward through both populations. For e = 0.4, 
nearly all pulsars including the 8.5 s PSR J2144-3933 (Young, Manchester & Johnston 1999) 
are able to produce pair cascades which could enable coherent radio emission. The lines 
decrease in slope as they move from longer to shorter periods and from small to large e 
as the particle acceleration becomes more limited by curvature radiation reaction. For the 
shortest period MSPs, the acceleration is completely radiation-reaction limited (Luo et al. 
2000) which causes a slower increase in electron Lorentz factor and impedes the development 
of pair cascades. The envelope of the observed pulsar population seems to be consistent with 
the shape predicted from this limit. 



4. Discussion 

We have found that distortions of a dipole magnetic field that produce small offsets of 
the PC can produce large increases in the multiplicity of pulsar pair cascades and allow pair 
cascades in pulsars not able to initiate cascades in centered dipole fields. Previous estimates 
of the effect of offset vacuum dipoles on pair death lines included only the decreased field 
line radius of curvature, concluding that dipole offsets of (0.7 - 0.8)i? ns (Arons 1997) or even 
0.95 R ns (Medin & Lai 2010), that are large fractions of a stellar radius, are required. Our 
calculation also takes into account the change in E\\ and the particle acceleration energy, 
which is by far the strongest effect. Since this effect is sensitive to effective offsets that 
are fractions of a PC radius, we show that many older pulsars and MSPs would be able to 
sustain robust pair cascades with dipole offsets that are small fractions of a stellar radius. 
Although distortions of the dipole field could be intrinsic to the NS through asymmetries 
in the interior currents, present either from birth or as a result of spin-up or spin-down 
evolution, magnetic dipole fields that are distorted by rotation or currents also produce 
offset PCs. The sweepback of field lines of a retarded vacuum dipole (Deutsch 1955) causes 
an offset of the PC towards the trailing side, opposite the direction of rotation (Dyks & 
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Harding 2004), and the force- free magnetosphere, with a similar sweepback of field lines, has 
even greater offset of the PC (Bai & Spitkovsky 2010). 

We will mention here a number of implications of our results. For high enough offsets, 
the entire radio MSP population could generate pair cascades. For the offset of e ~ 0.6 
implied by the apparent shift of ~ 1 km of the heated PC of PSR J0437-4715 (Bogdanov 
et al. 2007), with P = 5.75 ms and P = 1.4 x 10~ 20 , a pair cascade multiplicity of about 
10 2 is possible (cf. Table 2). The ability of MSPs with distorted dipoles to generate high- 
multiplicity pair cascades that can screen the E» over most of the open field could possibly 
account for the prevalence of sharply-peaked light curves of 7-ray MSPs that would require 
narrow acceleration gaps. Higher pair multiplicity in both normal and MSPs could also 
enhance their contribution to local Galactic cosmic-ray positrons. However, we find that 
the saturation in pair cascade multiplicity with increasing particle energy limits the extent 
to which a dipole distortion can produce even higher multiplicity. Therefore, even large 
values of e cannot produce the levels of multiplicity in young pulsars required to produce 
synchrotron emission observed from PWNe. We find for example that the Crab pulsar, with 
P = 0.033 s and B ~ 3 x 10 12 G could only produce a pair multiplicity up to ~ 3 x 10 4 
for e = 0.4 corresponding to an offset from Eqn (3) of 6% of the stellar radius. A pair 
multiplicity of ~ 10 6 in PSR J0737-3039A would also not be achievable with any degree of 
offset. 

Distorted dipole fields with offset PCs will generally introduce asymmetries in the pulsar 
emission as well as in the pair multiplicity. Since the particle acceleration and field line 
curvature vary over the PC, the pair multiplicity will be larger over one half of the PC in 
the direction of the offset. Pulsars that are below the death line for centered dipole fields, 
including many MSPs, will have pair cascades only on that side of the PC. The variation 
in pair multiplicity across the PC could result in asymmetric radio emission depending on 
how the coherent mechanism is related to the multiplicity. The higher accelerating field 
and voltage on the offset side of the PC will produce asymmetric heating of the PC and 
thus inhomogeneous thermal X-ray emission. The structure and energetics of the proposed 
slot gaps (SG) that form between the boundary of the open magnetic field and the upward 
curving PFF (Arons & Scharlemann 1979), and can accelerate particles to high altitude 
(Muslimov & Harding 2004), could be strongly affected by a distorted dipole. The particle 
Lorentz factor 7 in the SG , which is expected to reach curvature radiation-reaction limit 
such that 7 oc E^ 4 , will be larger on one side of the PC, producing CR emission power 
proportional to E\\ that is larger than for a dipole field. An azimuthal asymmetry of both 
the radiation power and width of the SG would change both the 7-ray luminosity and light 
curves. A further application of these results is in pulsar population studies. The decrease 
in the observed radio pulsar population just above the death line as pulsars spin-down and 
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evolve from left to right across the P-P diagram has been difficult to explain. But a range 
of dipole offsets across the population would produce a spread of pair/radio death lines and 
thus produce a death valley. It is interesting that the no-offset death line shown in Figure 1 
occurs near the densest part of the normal pulsar population. We will explore many of these 
consequences of increased pair multiplicity in distorted dipole fields in future studies. 
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J. Arons, A. Timokhin, O. DeJager and A. Spitkovsky, provided stimulation for this work. 
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Program, the Fermi Guest Investigator Program and the Universities Space Research As- 
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Table 1: Pair Cascade Multiplicity for Normal Pulsars 
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For x — 60°, <p — 270°, £ = 0.5 and neutron star parameters M = 1.45, R QS = 10 km and 
I = 1.13 x 10 45 gcm 2 . 
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Table 2: Pair Cascade Multiplicity for Millisecond Pulsars 
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For x — 45°, 4> = 270°, £ = 0.5 and neutron star parameters M = 2.15, R ns = 9.9 km and 
/ = 1.56 x 10 45 gcm 2 . 




Fig. 1. — Field lines of distorted magnetic dipole having an offset polar cap in the x-z plane 
and offset parameter e = 0.2. 
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Fig. 2. — Death lines in the P-P diagram for pair production by curvature radiation, for 
different values of the offset parameter e and inclination angle \ = 60°. The NS radius R QS , 
mass M ns (in Solar mass units) and moment of inertia J45 = J/10 45 gcm 2 refer to different 
NS equations of state used for normal and millisecond pulsars and are described in the 
text. Radio pulsars with measured P from the ATNF catalog (Manchester et al. 2005, 



http://www.atnf.csiro.au/research/pulsar/psrcat) are plotted as black dots. 



